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ABSTRACT 

We consider electrostatic plasma modes along the open field lines of a rotating 
neutron star. Goldreich-Julian charge density in general relativity is analyzed for the 
neutron star with zero inclination. It is found that the charge density is maximum at 
the polar cap and it remains almost same in certain extended region of the pole. For 
a steady state Goldreich-Julian charge density we found the usual plasma oscillation 
along the field lines; plasma frequency resembles to the gravitational redshift close to 
the Schwarzschild radius. We study the nonlinear plasma mode along the field lines. 
From the system of equations under general relativity, a second order differential 
equation is derived. The equation contains a term which describes the growing plasma 
modes near Schwarzschild radius in a black hole environment. The term vanishes with 
the distance far away from the gravitating object. For initially zero potential and field 
on the surface of a neutron star, Goldreich-Julian charge density is found to create the 
plasma mode, which is enhanced and propagates almost without damping along the 
open field lines. We briefly outline our plan to extend the work for studying soliton 
propagation along the open field lines of strongly gravitating objects. 

Subject headings: MHD - plasmas - pulsars: plasma: general - relativity - stars: 
neutron 



-3- 



1. INTRODUCTION 

Study of plasma modes in the neutron star or black hole environments is related with the 
investigation of radio emissions coming from these sources (see, e.g. Buzzi et al. 1995, Mofiz 1997 
and the references therein). Radio pulsars which are rotating neutron stars with spin periods 
ranging from ms to 5 s, are characterized by surface magnetic fields of the order of 10 12 G, radii of 
about 10 km and central densities in excess of 10 14 g ■ cm 3 , and so are purely gravitating objects. 
A spining magnetized neutron star generates huge potential differences between different parts of 
its surface (Goldreich k, Julian 1969 ). The cascade generation of electron-positron plasmas in the 
polar cap region (Sturrock 1971, Ruderman & Sutherland 1975) means that the magnetosphere 
of a neutron star is filled with plasma - screening the longitudinal electric field. This screening 
results in the corotation of plasma with a star. Such a rotation is not possible outside the light 
cylinder, thus it forms essentially different groups of field lines: closed i. e. those returning the 
stellar surface, and open, i.e. those crossing the light cylinder and going to infinity As a result, 
plasma may leave the neutron star along the open field lines. The charges along the field lines 
create plasma modes which may be related with the pulsar radiation and with its microstructures. 

Our study of plasma modes along the field lines is boosted by the pioneering works of 
Goldreich k Julian (1969), Sturrock (1971), Mestel (1971), Ruderman & Sutherland (1975) and 
Arons & Scharleman (1979). The subsequent achievements and some new ideas are reviewed 
by Arons (1991), Michel (1991), Mestel (1992) and Muslimov & Harding (1997). Although a 
self consistent pulsar magnetosphere theory is yet to developed, the analysis of plasma modes 
in the pulsar magnetosphere based on the above mentioned works provides firm grounds for the 
construction of such a model. 

In this paper, we attempt to extend Muslimov & Harding work (1997) to study plasma 
modes along the open field lines of a rotating neutron star. In §2 general relativistic equations 
describing the electrodynamics of a rotating neutron star are formulated. The equations are 
rewritten in the frame of reference corotating with the neutron star. We deduce the general system 
of equations governing the electrostatic modes in the pulsar magnetosphere. A detail analysis of 
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Goldreich-Julian charge density in general relativity in done in §3. It is shown that the charge 
density exponentially decays with the distance away from the surface of the star while it has a 
periodical dependence on the polar angle along the surface. The field is maximum at the polar cap 
region and it remains almost same in certain extended region in the pole. In §4 we study the linear 
plasma modes along the open field lines. A general equation governing electrostatic potential is 
derived. For a steady state Goldreich-Julian charge density, the usual plasma oscillation along the 
field lines is found. Plasma frequency resembles to the gravitational redshift close to Schwarzschild 
radius while at a large distance from the gravitational radius, it is the usual plasma oscillation 
along the field lines. In §5 we study the nonlinear plasma modes along the field lines. From the 
system equations under general relativity, a second order differential equation is derived. The 
equation contains a term which describes the growing plasma mode near the Schwarzschild radius 
of a neutron star or a black hole. The term vanishes with the distance far away of the gravitating 
object. The equation is solved numerically subjected to appropriate boundary conditions. It is 
found that Goldreich-Julian charge density creates the initial field on the surface of the star which 
is enhanced near the gravitational radius and almost without damping propagates along the open 
field lines. In §6 we conclude our findings and discuss them for further investigations. 

2. GENERAL RELATIVISTIC ELECTRODYNAMIC EQUATIONS IN THE 

COROTATING FRAME OF REFERENCE 

Recently Muslimov and Harding (1997) derived the general relativistic electrodynamic 
equations for a neutron star in the corotating frame of reference. It is noted that the effects of 
general relativity are very important : the dragging of inertial frames of reference significantly 
affects the electric field generated in the vicinity of a rotating magnetized neutron star, while the 
static part of the gravitational field results in additional enhancement of electric and magnetic 
fields near a star. 

The metric of an asymptotically flat, stationary, axially symmetric spacetime around a 
rotating gravitating body (see, e.g. Landau & Lifshitz 1975) is considered. In spherical polar 
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coordinates x° = ct, x 1 = r, x 2 = 9 and x 3 = <p, we have 

ds 2 = A 2 {cdt) 2 - B 2 {dr) 2 - C 2 {d6) 2 - D 2 {d(f) - codt) 2 , (1) 

where A = B~ l = (1 — r g jr) x l 2 is the gravitational redshift function, C = r,D = rsin6>, 
r g = 2GM/c 2 is the gravitational radius of body (neutron star) of mass M, J is the angular 
momentum of a neutron star, c is the speed of light and G is the gravitational constant. The 
metric in equation (1) is the approximation of Kerr metric when the ratio J/Mcr g is small. The 
presence of the nondiagonal component in metric in equation (1) results in the well known effect 
of dragging of inertial frames of reference (the Lense-Thirring effect) with the angular velocity 

2GJ 

- = ^3- (2) 

The metric in equation (1) can be transformed to the frame of reference corotating with a 
neutron star: 

ds 2 = A 2 (cdt) 2 - B 2 (dr) 2 - c 2 {d6) 2 - D 2 (d(j) - tldt) 2 , (3) 

by transformations t' = t,r = r, 6' = 6,<p' = (p — Slot. Here Q = to — £lo, where f^o is the angular 
velocity of rotation of star relative to the distant observer. 

The zero angular momentum observer (ZAMO; see, e.g. Thorne, Price & Macdonald 1986) 
has the four-velocity 

i n 



e u { , ,0,0, , }; ^{-^1-7^,0,0,0}. (4) 

yjl-rg/r cyjl-rg/r 

Then the general-relativistic Maxwell equations for observer (eq. [4] ) in the metric (eq. [3] ) take 
the form 

V B = , (5a) 



1 <9R 

Vx(aE-( W -fi )xB) = --- , (56) 
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V • E = 4vrp, (5c) 



1 <9E 47r 

V x (qB) + V x ((u - n„) x E) - (w - fi„) (V • E) = -— + — aj , (5d) 

c at c 



where (O x B) a = e^^epQ^, (u - ft ) a = {o,0,0, [(w - n o )"in0] / Cy Jl-r g /r}. 

Similarly, we may write the charge continuity equation in the above mentioned frame as 

dp ( k 



. )t + ( — -!) n m- Vp + V ■ (oj) = 0. (6) 



Finally, the equation of motion of a charged particle is 



1 dp ( v , 

—jr = mjg + q E + - x B +f . (7) 



Here a = (1 — r^/r) 1 / 2 (= A, as denoted in the metric in equation (1) ), the parameter 
k = 2r g R 2 /5, R is the radius a neutron star, p = J2 S n sQs,3 = J2 S n sQsV s , v s is the velocity, q s is 
the charge of particle, n s is the particle number density, and summation is over all the species; 
p = 7717V is the momentum of the particle, 7 = (1 — v 2 /c 2 ) -1 / 2 is the Lorentz factor, m is the rest 
mass of the particle, f is an external force other than electromagnetic, and g is the gravitational 
acceleration. All electrodynamic quantities as magnetic field B, electric field E, conduction current 
j, and charge density p in these equations are such as measured by ZAMO (eq.[4]). Gradient, curl 
and divergence are taken along the curvilinear coordinate 

a 3 lid 1 Id 

ef = Ae r = A—, e§ = -e e = -—, e ; = — — e = — ^^7, (1) 
or ° r r 00 9 r sin r sin aq> 

m = r sin Oe^ is the Killing vector, responsible for axial symmetry. 

Assuming ^ = in equation (5b) (i.e. considering magnetic field of a neutron star is 
stationary in the corotating frame), from equation (5b) we get 

qE-(w-O )xB = -V*, (8) 
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where $ is a scalar electrostatic potential. 



Taking the divergence of equation (8) and making use of equation (5c), we get 

V • ( -V* + -(O - w) x = -4vrp . 
\a a J 

Equation (9) can be written as (see Muslimov & Tsygan 1986; Beskin 1990) 



(9) 



V- (-V*) = -4^(p-p GJ ) 



a 



where 



PGJ = -i-V • (i(no - «) x B) = -^V • - £)Oo x B 

is the relativistic analog of the Goldreich-Julian (1969) charge density. 



Finally, the equation of motion of the charged particle is 



d I h 
0t + 



1 J n ■ V + av • V 



-gV$, 



(10) 

(11) 



(12) 



where the gravitational acceleration g and the nonelectromagnetic force f are justifiably ignored. 



3. GOLDREICH-JULIAN CHARGE DENSITY IN GENERAL RELATIVITY 

In a pioneering work, Goldreich & Julian (1969) have shown that a strongly magnetized, highly 
conducting neutron star, rotating about the magnetic axis, would spontaneously build up a charged 
magnetosphere. The essence of the argument is that it imposes a charge magnetosphere which 
are subject to enormous unbalanced electric forces parallel to the magnetic field. Goldreich and 
Julian hypothesized that a far better approximation for the magnetosphere would be shorting-out 
of the component of E along B by charges originating in the star. The magnetospheric charges 
that maintain E • B = are themselves subject to the E x B drift which sets them into corotation 
with the star. Here, we analyze Goldreich-Julian charge density in general relativity. 

Assuming zero inclination of the rotating star with the magnetic axis we consider 
B = {B r ,Bg,0}. The components B r ,Bg in this case were first derived by Ginzburg & Ozernoy 
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(1964). Later on, similar expressions were derived in a number of papers (see, e.g. Wasserman & 
Shapiro 1983; Muslimov& Tsygan 1986): 



B r — 



2cos# 



(13a) 



where 



Ba 



sin6 



(136) 



ff \ 3r3 



s L 



In|l-^+- + U- 

r / r 2 \ r 



(14) 



V>(r) 



3r 2 



1 - ^ r 



( + 2 - In ( 1 - ^ ] + 1 



r 



and /u is the magnetic dipole moment of a neutron star. 



(15) 



We perform a detail calculation of pcj with the magnetic field of a rotating neutron star 
given by equations (13a) and (13b), respectively. From equation (11), we find 



PG.J 



iircr 2 sin 6 



k \ r 3 sin 2 9Ba 



1 — r g /r 



1 - s 



r 2 sin 2 6B r 



r 3 J i _ r 9 /r 



The calculation shows that 



(16) 



with 



p GJ (r,#) = fi ? i(r)sin 2 ^-F 2 (f)(sin 2 ^-2cos 2 0) 
F 1 (r-)=r-{((l-|){ ? | T - ? i p + 21„(l-i)} 



(17) 



(18) 
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Tg. Asymptotically as r/r g — > oo functions F] 



(19) 



Thus Goldreich-Julian space charge has two purely general relativistic contributions, one 
is due to the Schwarzschild gravitoelectric parameter r g and the second one is due to the 
gravitomagnetic Kerr parameter f3. They have different dependence on r as 1/r and 1/r 3 , 
respectively. It is meant that near the surface of the star the gravitomagnetic term is in concurent 
with the gravitoelectric one. But in the distance from the surface of the star which is comparable 
with its radius R the gravitomagnetic term is ignorable small. 

We plot Fi(f) and ^(f) for (3 = 0.1. The dependence of these functions on f is shown in 
Fig.fl] and Fig [2], respectively. The Goldreich-Julian charge density under general relativity is 
shown in Fig. [3]. 

By least-squares fitting of the curve at 6 = 0, we find that Goldreich-Julian charge density 
decays with the distance away rom the star as follows: 



10.5053 



5.05692 



1.06093 



- 0.084179 . 



(20) 



PGJ = 



.3 



+ 



r 



r 



The charge density is maximum at the polar cap region and it remains almost same in certatin 
extented region in the pole. It is to be noted that the expression for pcj obtained by Muslimov 
and Harding (1997) shows the similar results. 
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Fig.[l]. Goldreich- Julian charge density in general relativity; F\ as function of r. 
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Fig. [2]. Goldreich- Julian charge density in general relativity; as a function of r 
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Fig. [3]. Goldrech- Julian charge density in general relativity pc,j as function of r. Thick line 
corresponds to 6 = 0, thin line to 8 = tt/4, and broken line to 8 = tt/3. 

4. LINEAR PLASMA MODES ALONG THE OPEN FIELD LINES 

The theory of cascade generation of electron-positron plasma at the polar cap region of a 
rotating plasma is developed by Ruderman & Sutherland (1975). According to the theoretical 
model, due to escape of charge particles along the open field lines , a polar potential gap is 
produced which continuously breaks down by forming electron positron pair on a time scale of 
a few microseconds. A photon of energy greater than 2mc 2 produces an electron-positron pair. 
The electric field of the gap accelerates the positron out of the gap and accelerates the electron 
towards the stellar surface. The electron moves along a curved magnetic field line and radiate an 
energetic photon which goes on to produce a pair as it has a sufficient component of momentum 
perpendicular to the magnetic field. Recently, Zhy & Ruderman (1997) explained the e — e + pair 
production from a Crab-like pulsar. Electrons and positrons are accelerated in opposite directions 
to extremely high energies. The Lorentz factor 7 of the primary electron and positron is given by 
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where This cascade of pair production, acceleration of electrons and positrons along curved field 
lines, curvature radiation-pair production results in a "spark" break down of the gap. 

Assuming a steady state thermodynamically equilibrium plasma state in the polar cap region, 
we study the linear plasma modes along the field lines. From the system of equations (10), (6) and 
(12), we derive the following linearized equations: 



V • ( 1 Vd>) = -47r(]T q s 8n s - p GJ ) (20) 



) + V ■ U/T^v.) = 0, (21) 
at' no V r 



£ = < 22 > 

of m 

where d/dt' = d/dt + - l) Slo m • V is the global time derivative along ZAMO trajectories, 
s(= e, e + ) is the plasma species, 5n s is the density fluctuation of the plasma species and no is 
the equilibrium plasma density and pcj is the Goldreich-Julian charge density as defined by the 
equation (17). The system of equations (20)-(22) is equivalent to the following equation: 

% (v • (-/=^ V<I> ) " 4 ^)j + V • {u P0 ^Jl^V^) = 0, 



(23) 



dt' 

where u>p = 8imQe 2 /m. 

Now, by defining the electric field arising from charge separation and the corotational electric 



field, which is the source of pcj 



E = J— V$, pcj = — V • E c , (24) 



from equation (23), we find 



dt 



a -[V-(E + E c )]+V 



= 0, (25) 



-13- 



which gives 







-E + ^Jl- 



E 



For |pf = we may write the solution of equation (26) as 



(26) 



E = Eq exp ^ —iu> po \j 1 -t' 



(26) 



From the above solution, we find that the plasma frequency in general relativity now is 
defined as 

"J = < i 1 ~ 7) ( 28 ) 
which is equivalent to the gravitational redshift of the oscillation. Fig. [4 ] shows the dependence 
of plasma oscillation on the distance away from the gravitational radius of the star. 




Fig. [4]. Plasma frequency uj p (r) as a function of r; gravitational redshift near the Schwarzschild 
radius. 

The global time derivative along ZAMO trajectories is defined as 

w ~ m + b " n "W (29> 



Thus, we may define 



(30) 
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and hence the solution of linear plasma mode is 



E(t, r, 0) = E exp { -iu Po x /l - -^-{t + — _ ^ 



(31) 



Introducing the dimensionless quantities e = J^,t = w po t,x = j- g ,X = t^->^ = 7% from the 
equation (31), we find 

' " ' ' ' (32). 



e(r, x, 4>) = sin 



1 




Fig. [5]. Linear plasma mode in general relativity; electrostatic field e(r, 4>) as a function of r. 

We do some analysis of the linear electrostatic modes around a rotating neutron star. First, 
we consider the Schwarzschild radius equal to half the radius of the neutron star, i.e., r g = R/2. 
Then we find 5 = 8/5. Considering relatively dense plasma, we put \ = 10 8 - We consider a fixed 
azimuthal angle (f) = it/ 4. For r = we plot the field e(f), which is shown in Figure 5. We find 
that the electrostatic field generated by Goldreich-Julian charge density is maximum near the star 
surface and falls quickly from the star. 
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5. NONLINEAR PLASMA MODES ALONG THE FIELD LINES 

Now, we consider nonlinear plasma modes along the open field lines around a rotating neutron 
star. The system of equations governing the nonlinear modes can be written as 



d_ 

di' 



7 (n s )+V ■(^l- r -^n s v s )=0, (33) 



757 + \ /l-^v.-V v s = -fVd>, (34) 



V • ( j± V$) = -4vr(^ n s q s - p GJ ) . (35) 

For the simplicity, we consider v s .V « w sr • Vd/dr (i.e. one dimensional wave propagation 
along r) and introduce a moving frame rj = r — Vt', where V is a constant. In the considered 
moving frame, from equations (33) and (34), we get 



Using equations (36) and (37), in equation (35), we derive the nonlinear equation for the 
plasma mode along the field line of the rotating neutron star: 

<?<$> r g d<$> , ujUl - r f) $ 



+ 



drj 2 2r) 2 {1 _r i)drj y 2 1 _ 2(1 _^ )( ^ )2+ 1 (1 _ ^2(^)4 



= 4tt Jl - ^p G j . (38) 



Now introducing dimensionless quantities 
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$ = 



LO, 



mv 



Pa 



V 



(39) 



we write the equation (38) in dimensionless form: 



+ ^o(l-^)- 



d?<f> 1 

~drf ~ 277(77 - 1) dfj ' ~ po ^ r7 ; i_2(l-I)$ 2 + i(i_I)2$4 



where 



F — — 

2F 2 



■ * 1- 



77/ 



Fi (77) sin 2 9 - F 2 (77) (sin 2 9-2 cos 2 0) 



(40) 



(41) 



Here Q. c = eBo/mc; Fi(rj) and ^2(77) are determined by the equations (18) and (19), respectively. 




Fig. [6]. Nonlinear plasma mode in general relativity; propagation of plasma oscillation $(77) 
near the surface of a neutron star; numerical solution of the equation (40) with boundary condition 
$(2) = 0, $'(2) = 0. 

We numerically solve the equation (40) in the polar cap region (9 s=s 0) of a neutron star 
subject to the appropriate boundary conditions. Following Goldreich & Julian (1969) and 
Muslimov & Harding (1997), we assume that the surface of a polar cap and that formed by the 
last open field lines can be treated as electric equipotentials. We therefore adopt the condition 
$(r = R) = 0. Second, we require that the steady state component of electric field parallel to 
magnetic field vanishes at the polar cap surface, i.e. d&(r = R)/drj = . By considering r g = R/2 
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for a neutron star, we write the boundary conditions as: $(2) = and $'(2) = 0. The solution of 
the equation (40) with the mentioned boundary condition is shown graphically in Fig. [6]. We find 
that Goldrech-Julian charge density creates the initial potential on the surface. Near the radius 
the potential is enhanced and it propagates almost without damping along the field lines. 

6. DISCUSSION AND CONCLUSION 

We study the electrostatic plasma modes along the open field lines of a rotating neutron star. 
The dragging of inertial frame and the effect of general relativity is fully considered in this study. 
We perform a detailed analysis of Goldreich-Julian charge density in general relativity. Since 
pulsars having smaller obliquity have larger accelerating drops and this favored for 7 -ray pulsar 
emissions (Muslimov 1995) and it supports the single pole 7 - ray pulsar models (Daugherty & 
Harding 1994, 1996; Dermer & Sterner 1994), we confine our analysis in the zero inclination of the 
rotating neutron star. As pulsar radiation takes place in the plasma environment or the radiation 
passes through a plasma media , we consider the electrostatic plasma modes along the open field 
lines. We study both the linear and nonlinear modes in the neutron star or black hole plasma 
environment. Our general conclusion from the above analysis may be summarized as follows: 

1. Goldreich-Julian charge density is maximum in the polar cap region and remains almost 
same in a certain extended region of the pole. The charge density exponentially decays with the 
distance away from the surface of the star. 

2. Plasma oscillation along the field lines resembles to the gravitational redshift near 
Schwarzschild radius. 

3. Plasma modes grows near the gravitational radius in the black hole environment. 

4. For initially zero field on the surface of a rotating neutron star, Goldreich-Julian charge 
density, which is enhanced near the surface and propagates almost without damping along the 
open field lines, creates the plasma modes. 
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For further study of plasma dynamics in the neutron star or black environment we plan to 
extend our earlier investigations on solitons (see, Mofiz 1989, 1990, 1993 & Mofiz et al. 1985, 
1995) propagation along the open field lines of strongly gravitating objects. 
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